Negativity as the Entanglement Measure to Probe the Kondo Regime in the Spin-Chain Kondo 

Model 
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We study the entanglement of an impurity at one end of a spin chain with a block of spins using negativity as 
a true measure of entanglement to characterize the unique features of the gapless Kondo regime in the spin chain 
Kondo model. For this spin chain in the Kondo regime we determine- with a true entanglement measure- the 
spatial extent of the Kondo screening cloud, we propose an ansatz for its ground state and demonstrate that the 
impurity spin is indeed maximally entangled with the cloud. To better evidence the peculiarities of the Kondo 
regime, we carry a parallel analysis of the entanglement properties of the Kondo spin chain model in the gapped 
dimerised regime. Our study shows how a genuine entanglement measure stemming from quantum information 
theory can fully characterize also non perturbative regimes accessible to certain condensed matter systems. 



I. INTRODUCTION 

The investigation of entanglement or the truly "quantum" 
correlations inherent in many -body condensed matter systems 
is currently a topic of intense activit)^^. This emerging area 
aims at characterizing many-body states using tools and mea- 
sures developed in quantum information. Till date, most in- 
vestigations have focussed on either the entanglement be- 
tween individual elements, such as single spins, or the entan- 
glement between two complementary blocks of a many-body 
system. The former entanglement is generically non-zero only 
between nearest or next to nearest neighbors^--. For comple- 
mentary blocks, the whole system is in a pure state and the 
von Neumann entropy is a permissible measure of the entan- 
glement. In this context much interest has been evoked by 
conventional gapless phases, where, due to the absence of an 
intrinsic length scale, the von Neumann entropy diverges with 
the size of the blocks^^. In this backdrop, it is timely to in- 
vestigate the entanglement in gapless regimes of a many body 
system for which its true form and amount may not be char- 
acterized through the entanglement between individual spins 
or complementary blocks. Kondo systems are ideal candidate 
for these investigations^iS. 

Kondo systemsiili are expected to be very distinctive in the 
context of entanglement for at least two reasons: (a) Despite 
being "gapless", they support the emergence of a length scale 
^- the so called Kondo screening lengthiiS- which should be 
reflected in the entanglement, making it markedly different 
from that in the more conventional gapless models studied so 
far; (b) They are expected to have a more exotic /on7z of entan- 
glement than the widely studied spin-spin and complementary 
block entanglements. Indeed, in Kondo systems, the impurity 
spin is expected to be mostly entangled with only a specific 
block of the whole system. This is, of course, merely an intu- 
ition which needs to be quantitatively verified with a genuine 
measure of entanglement: this is indeed the task accomplished 
in this paper, where we provide the only characterization of 
the Kondo regime based entirely on a true measure of entan- 
glement. 

The simplest Kondo mode l'"'" describes a single impurity 
spin interacting with the conduction electrons in a metal; the 



ground state is a highly nontrivial many body state in which 
the impurity spin is screened by conduction electrons in a 
large orbital of size ^. Many physical observables vary on 
the characteristic length scale ^, which is a well defined func- 
tion of the Kondo coupling '*'. Determining the spatial extent 
of the Kondo cloud has been so far a challenging problem re- 
peatedly addressed by various mean a^-'^-'^ . This includes an 
investigation which introduces a quantity called "impurity en- 
tanglement entropy" which, however, is not a measure of en- 
tanglement and cannot quantify the entanglement within the 
systeniiS. Recentlj^ii, it has been pointed out that the uni- 
versal low energy long distance physics of this Kondo model 
arises also in a spin chain when a magnetic impurity is coupled 
to the end of a gapless Heisenberg anti-ferromagnetic Ji — J2 
spin 1/2 chain, where Ji ( J2) is the (next) nearest neighbor 
coupling. When J 2 exceeds a critical value, the spin chain en- 
ters a gapped dimerized regime and its relation to the Kondo 
model breaks down. Namely, for < J2 < J| = 0.2412, the 
spin system is gapless and it supports a Kondo regimai^. For 
J2 > the system enters the gapped dimer regime, where 
the ground state takes a dimerised form; at the Majumdar- 
Ghosbi^ point ( J2 = 0.5), the ground state becomes just a 
tensor product of singlets. For J2 > 0.5, incommensurability 
effects'- emerge. 

Our aim in this paper is to use a true measure of entan- 
glement to fully characterize the unique features of the gap- 
less Kondo regime in the spin chain Kondo model. Namely, 
for this spin chain in the Kondo regime: (i) we demonstrate 
that the impurity spin is indeed maximally entangled with 
the Kondo cloud; (ii) we determine the spatial extent of the 
Kondo screening length ^ using only an entanglement mea- 
sure; (iii) we motivate an ansatz for its ground state in the 
Kondo regime; (iv) we evidence the scaling of a true measure 
of entanglement as pertinent parameters are varied. In order 
to accomplish these tasks we device a DMRG approach en- 
abling to investigate the entanglement between a single spin 
and a pertinent block of the chain, which may be applied in 
other contexts. Finally, to better evidence the unique prop- 
erties of the entanglement in the Kondo regime we carry a 
parallel analysis of the entanglement properties of this model 
in the gapped dimerised regime. Using a true measure of en- 
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FIG. 1 : (Color online) a) Kondo Spin chain with next nearest neigh- 
bor Heisenberg interaction with one impurity at one end. b) The 
chain is divided into three parts, an impurity, a block A and a block 
B. Entanglement is computed between the impurity and block B. 

tanglement to determine ^ enables to exploit the peculiarities 
of the Kondo regime of a spin chain to generate long range 
distance independent entanglement usable for quantum com- 
munication tasksii. 

A true measure of entanglement should satisfy a set of pos- 
tulates - for example, it should be non-increasing under lo- 
cal actions: such a genuine measure does exist for two sub- 
systems of arbitrary size even when their combined state is 
mixed, as it happens in Kondo systems. This measure is the 
negativity^ and it has been successfully used to quantify the 
entanglement in a harmonic chai n and between distant re- 
gions of critical systems^i^. For bipartite systems, negativity 
is defined as E = \ai\ — 1, where denote the eigen- 
values of the partial transpose of the whole density matrix of 
the system with respect to one of the two subsets of the given 
partition and |...| is the absolute valuai^. 

The paper is organized as follows: In section (|II]i we de- 
fine an entanglement healing length for the spin chain Kondo 
model; section ( illli ) explains the DMRG-based approach we 
devised in order to compute its entanglement properties. In 
section dlVI ) we show the remarkable scaling of a true measure 
of entanglement (i.e. negativity) in the Kondo regime attain- 
able by the Kondo spin chain when < J2 < J2 = 0.2412; 
in addition, we motivate an ansatz for the ground state of this 
chain in the Kondo regime. Section (|V]) is devoted to a sum- 
mary of our results and to a few concluding remarks. 

II. MEASURING THE ENTANGLEMENT HEALING 
LENGTH 

The spin chain Kondo modelii is defined by the Hamilto- 
nian 

N-l N-2 
H — j'(tTi.tT2 + J^2CTl.CT3) + ^ CTi.fJi+i + J2 ^ Cri.tTi+2, 

i=2 i=2 

(1) 



FIG. 2: (Color online) a) DMRG representation of the state of the 
chain keeps two intermediate spins in ordinary computational basis 
and the left and right blocks in a truncated DMRG basis, b) The in- 
termediate spin next to the impurity is traced out from the density 
matrix of the chain. This tracing is equivalent to adding the traced 
out spin to block A. c) The basis of the right block of DMRG rep- 
resentation is transformed so that a single spin in the left side of the 
right block is represented in the computational basis while the state 
of the new right block is given in a DMRG basis. 



where ai = (erf , ct^, erf) is a vector of Pauli operators at site 
i, N is the total length of the chain, J2 is the next nearest 
neighbor coupling and the nearest neighbor coupling Ji has 
been normalized to 1. The impurity spin, located at one end 
of the chain, is accounted for by weaker couplings to the rest 
of the system; in the following, see Fig. [T^, both couplings Ji 
and J2 are weakened by the same factor J', which quantifies 
then the impurity strength. 

To study the entanglement of the ground state we divide- 
see Fig. [it) - all the spins of the chain in three different groups: 
the impurity spin, block A, which contains the L spins next to 
the impurity {L = 0,1, N — 1) and block B formed by the 
remaining N — L — 1 spins. We use negativity to fully char- 
acterize the entanglement between the impurity and block B 
in both the gapless Kondo and the gapped dimerised regimes. 
We determine the size of the block A when the entanglement 
between the impurity and block B is almost zero; by this pro- 
cedure we measure an Entanglement Healing Length (EHL) 
L*, i.e. the length of the block A which is maximally entan- 
gled with the impurity. We show that, in the gapless Kondo 
regime, EHL scales with the strength of the impurity coupling 
just as the Kondo screening length, ^, does. Thus, in the gap- 
less regime of the Kondo spin chain, our approach yields a 
method to detect the Kondo screening length^iSiil based on a 
true measure of entanglement. In addition, we find that entan- 
glement, as quantified by negativity, is a homogeneous func- 
tion of two ratios: N/L* and L/N, where L is the size of the 
block A, i.e. the block adjacent to the impurity, and N is the 
length of the whole chain. As a result, the entanglement in the 
Kondo regime is essentially unchanged if one rescales all the 
length scales with the EHL L*. 
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III. DMRG ANALYSIS OF ENTANGLEMENT IN THE 
SPIN CHAIN KONDO MODEL 



We use the Density Matrix Renormalization Group 
(DMRG^ approach to compute the ground state of the spin 
chain Kondo model. We analyze large chains, N = 250, to 
avoid finite size effects and take N to be even to avoid prob- 
lems arising from accidental degeneracies. In a DMRG ap- 
proach the ground state of the system is partitioned in terms 
of states of a left block, a right block (not to be confused with 
blocks A and B) and two intermediate spins as shown in Fig. 
|2^. The states of the intermediate spins are given in the com- 
putational (I t); I i)) basis, while the states of the both blocks 
are usually in some non-trivial truncated DMRG basis. In this 
approach one has several representations for the ground state 
which vary due to the number of spins in the left (right) block 
and it is possible to go from one representation to the other by 
applying pertinent operators on each block. The main issue 
in the DMRG is that the dimension of the left (right) block is 
kept constant independent of whatever spins are there in that 
block. To have a fixed dimension for the left (right) block 
we truncate the Hilbert space such that the amount of entan- 
glement between the two parts of the chain remains almost 
unchanged22,. To have a precise results we need to sweep all 
representations of the ground state for several times to get the 
proper basis for the left and the right blocks of all representa- 
tions. After some sweeps, when the ground state energy con- 
verges (we keep states for which the error of the energy would 
be less than 10^®), we pause to compute the entanglement. 
We take a representation of the ground state in which the left 
block contains just the single impurity spin and the right block 
contains iV — 3 spins: as a result, the single impurity spin is 
given in the computational basis and this allows us to com- 
pute the negativity later. From this DMRG state, we trace out 
the spins belonging to block A before computing the entan- 
glement between the impurity and block B since it is most 
convenient to compute the entanglement between the impu- 
rity and the block B: due to the entanglement monogamy, this 
provides an equivalent information about the entanglement of 
the impurity with the block A. Our tracing procedure starts 
with the density matrix of the ground state of the system in 
the representation shown in Fig. |2^; at this stage, the number 
of spins in the block A is zero (no spin has been traced out), all 
spins except the impurity belong to the block B, and the en- 
tanglement between the impurity and the block B is maximal 
(i.e. E = 1). Then, we trace out the intermediate spin next to 
the impurity as shown in Fig. [2}5; this amounts to putting that 
spin into block A. Finally, as shown in Fig. |2};, we transform 
the DMRG basis of the right block so as to put a single spin 
at the left of the right block in the computational basis, while 
the state of the new right block is given in a DMRG basis. As 
a consequence, the resulting density matrix has the exact form 
of Fig. and we can continue the procedure to trace one spin 
at each step (i.e., put more spins in the block A) and compute 
the entanglement between impurity and block B. 
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FIG. 3: (Color online) a) L* vs. 1/V J' for both Kondo (J2 = 0) 
and dimer regime ( J2 — 0.42). b) Entanglement vs. L/N for fixed 
N/L* — 4 when J2 = 0. c) Entanglement vs. L/N for fixed 
N/L* = 4 at the critical point J2 — J2. d) Entanglement vs. L/N 
for fixed N/L* = 4 in the dimer regime (J2 = 0.42). 



IV. SCALING OF NEGATIVITY AND ANSATZ FOR THE 
GROUND STATE IN THE KONDO REGIME 

We find that there is an EHL L* so that, for L > L*, the en- 
tanglement between the impurity and block B is almost zero: 
L* provides us with an estimate of the distance for which the 
impurity is mostly entangled with the spins contained in block 
A. For large chains (N > 200) in the Kondo regime, one 
finds that L* is almost independent of N and depends only on 
J'. In the Kondo regime, i.e. for J2 < J|, L* depends on 
J' just as the Kondo screening length ^ does^^; for small J', 
L* oc e"/^^, where a is a constant. We plot L* as a function 
of I/VJ' in Fig. [3^. In a semilogarithmic scale, the straight 
line plot exhibited in the Kondo regime ( J2 = 0) shows that 
L* may be regarded as the Kondo screening length. More- 
over, the nonlinearity of the same plot in the dimer regime 
(J2 = 0.42), especially for small J', shows that, here, no ex- 
ponential dependence on X/^/lP holds. 

We observe also a remarkable scaling of negativity in the 
Kondo regime. This scaling may be regarded as yet another 
independent evidence of the fact that L* is indeed the Kondo 
length ^. In general, the entanglement E between the impurity 
and block B is a function of the three independent variables, 
,/', L and iV which, due to the one to one correspondence be- 
tween J' and L*, can be written as E{L* ,L,N). We find 
that, in the Kondo regime, E = E{N/L*,L/N). To illus- 
trate this, we fix the ratio N/L* and plot the entanglement in 
terms of L/N for different values of J' (or equivalently L*) 
for J2 = (Fig.[3}5) and for J2 = J| (Fig. [3};). The complete 
coincidence of the two plots in Figs. [3J3 and c shows that, in 
the Kondo regime, the spin chain can be scaled in size with- 
out essentially affecting the entanglement as long as L* is also 
scaled. In the dimer regime the entanglement stays a function 
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of three independent variables, i.e. E = E{L* , L, N), and, as 
shown in Fig. |3}l, the entanglement does not scale with L*. In 
our approach, the Entanglement Healing Length L* may be 
evaluated in both the Kondo and the dimer regime: the scal- 
ing behavior, as well as the dependence of L* on J', discrimi- 
nates then between the very different entanglement properties 
exhibited by the spin chain Kondo model as J2 crosses J|. 

We defined L* such that there is no entanglement between 
the impurity and block B when block A is made of L* 
spins. Conventional wisdom based on previous renormaliza- 
tion group analysis suggests that, in both regimes, the impu- 
rity and the block A of length L* form a pure entangled state, 
while block B is also in a pure state. This is indeed approxi- 
mately true in the dimer regime (exactly true for J2 = 0.5) but 
it turns out to be dramatically different in the Kondo regime. 
To check this, we computed the von Neumann entropy of the 
block B when block A has L* spins and found it to be non 
zero. Thus, the blocks A and B are necessarily entangled in 
the Kondo regime as there is no entanglement between the 
impurity and B. In fact, after a distance L*, the impurity is 
"screened" i.e, the block B feels as if it is part of a conven- 
tional gapless chain and has a diverging von Neumann en- 
tropy. The Kondo cloud is maximally entangled with the im- 
purity as well as being significantly entangled with block B. 
Based on the above, a simple ansatz for the ground state \GS) 
in the Kondo regime is provided by 



|G5)=E' 



t)i4(J0)-u}|Lt(J0) 

V2 



^\R,{J')), (2) 



where a, are constants, {|iJ(J')>, iLjiJ'))} and {\R,{J'))} 
are sets of orthogonal states on the cloud and the remaining 
system, respectively. At the fixed point J' — > all spins 
except the impurity are included in \lJ{J')) and \Lf{J')). 
At J' — >■ 1, very few spins are contained in |lJ(J')) and 
\Lf{J')) while {\Ri{J'))} represents most of the chain. 

For what concerns the mere evaluation of the amount of en- 
tanglement as J' is varied, we first plot, in Fig. |4^, the negativ- 
ity as a function of L near the two fixed points, i.e. J' — > and 
J' — > 1, accessible in the Kondo regime: as expected, near 
J' — > (i.e, for large values of the Kondo screening length), 
the entanglement remains large for rather large values of L 
while, for J' ^ 1 (i.e. for a very small cloud) it decreases 
rapidly with L. Fig. (semilogarithmic) shows that, also 
at the extreme limits J' ^ and J' ^ 1, the entanglement 
decays exponentially with L since this a characteristic mark 
of the entanglement in the Kondo regime. This exponential 
decay of entanglement is absent in the dimer regime: Fig. |4j3 
shows that, here, only for rather large J', the entanglement 
decays exponentially with L while, for small J', the entangle- 
ment between the impurity and block B decays slower than 
an exponential as a function of L exhibiting even a plateau 
at short distances. This latter feature is evidenced in Fig. |4};, 
and is consistent with the emergence, for small J', of long 
range valence bonds between the impurity and far spins as a 
consequence of the onset of the dimerised ground stated. In 
fact, when J' is small, J2 J' becomes much less than J| and 
the impurity forms valence bonds with distant spins, while the 
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FIG. 4: (Color online) a) Entanglement vs. L for the two fixed points 
J' — ;> and J' — >■ 1 in the Kondo regime J2 = 0. b) Exponential 
decay of entanglement in terms of L in a chain of length A'^ = 250 
for J' = 0.6. c) Non-exponential decay for small J' in the dimer 
regime J2 = 0.42. 



Other spins, since for them J2 > J2, form a valence bond 
with their nearest neighbor to preserve the dimerised nature 
of the ground state: as a result, the impurity shares less entan- 
glement with nearby spins and fulfills its capacity of entangle- 
ment forming valence bonds with the more distant spins in the 
chain. 



V. SUMMARY AND CONCLUDING REMARKS 

To summarize, we analyzed the Kondo spin chain model 
from the viewpoint of a genuine entanglement measure, 
namely the negativity. This readily shows that the impurity 
spin is indeed maximally entangled with the Kondo cloud. 
We used negativity to provide an independent method to de- 
termine the Kondo screening length and to provide a charac- 
terization of the ground state of the Kondo spin chain in the 
Kondo regime. We found that, not only is the Kondo regime 
of this model distinct from the gapless phases probed to date 
using the von Neumann entropy, but the form of the entan- 
glement - a spin and a block in a mixed state - is also dis- 
tinctive. We deviced a DMRG approach enabling to compute 
the entanglement between the impurity and a block of spins 
located at the other side of the chain for different lengths of 
the block. We showed that, in the Kondo regime, the EHL L* 
scales with the impurity coupling J' just as the Kondo length 
does: in other words, the impurity, though not entangled with 
any individual spin, is shown to be entangled with the totality 
of the spins within the Kondo cloud- whose size is measured 
by L*- and disentangled from the rest. Our measure of the en- 
tanglement in the Kondo regime led us to formulate an ansatz 
for the ground state of the Kondo spin chain for J2 < -/I- 
Our approach also shows that, in the Kondo regime, the en- 
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tanglement scales exponentially with L/L* and that, in the 
gapped dimer regime, though it is still possible to define an 
EHL, the impurity-block entanglement is usually smaller and 
has no characteristic length scale. 
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